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I.  INTRODUCTION 


The  study  of  the  3-D  perturbed  motion  of  a  rotating  fluid  which  fills  a 
spinning  and  coning  cylinder  shows  that  the  solution  depends  on  the  cylinder 

aspect  ratio  A  =  c/a  and  the  Reynolds  number  Re  =  a2$/v,  where  c  =  half 

height,  a  =  radius,  $  =  spin  of  the  cylinder  and  v  =  kinematic  viscosity  of 
the  fluid.  The  forced  oscillation  problem  was  solved  for  the  pressure  in 
Reference  1  and  for  the  moment  exerted  by  the  fluid  on  the  cylinder  in  Refer¬ 
ence  2  assuming  linear,  viscous  perturbations  of  the  steady  state,  i.e.,  solid 

body  rotation.  A  modal  analysis,  matched  asymptotic  expansions  for  the  flow 
near  the  endwalls  and  an  expansion  in  outer  flow  spatial  eigenfunctions  were 
empl oyed. 

The  ranges  of  interest  of  the  parameters  are  0.5  <  A  <  5  and  1  <  Re 
<  106,  determined  by  the  values  for  which  experimental  data,  of  various 
sorts,  are  available.  The  theory  of  References  1  and  2  is  asymptotic  for 
Re  “  and  requires  the  solution  of  an  eigenvalue  (e.v.)  problem.  The 

complex  e.v.  are  denoted  by  X  and  the  eigenfunctions  are  sin  Xx,  cos  Xx.  The 

"basic  set,"  X^,  satisfies  X^  kir/2A  for  Re  ",  where  k  =  1,3,5...  Only 

this  set  was  considered  in  References  1-4. 

The  existence  of  another  eigenvalue,  X$,  where  X$  =  0(Rely/2)  and 

therefore  is  not  in  the  basic  set,  is  shown  here;  in  this  limit  the  coning 
frequency,  t,  is  fixed  and  r  *  +  1,3.  The  limits  of  all  the  e.v.'s  as  coning 

frequency  approaches  spin  frequency,  r  *  1 ,  with  Re  fixed,  are  also 

considered. 


1.  Gerber,  N.,  Sedney,  R.,  and  Bart  os ,  J.  M.,  "Pressure  Moment  on  a  Liquid- 
Filled  Projectile:  Solid  Body  Rotation,"  U,  S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland ,  ARBRL-TR-02422 ,  October 

1982.  (AD  A12GS67 ) 

2.  Gerber,  N.,  and  Sedney,  R.,  "Moment  on  a  Liquid-Filled  Spinning  and 
Nutating  Projectile:  Solid  Body  Rotation,"  U.  S.  Army  Ballistic  Research 
Laboratory ,  Aberdeen  Proving  Ground,  Maryland,  ARBRL-TR-0247  0 ,  February 

1983.  (AD  A125332) 

3.  Jedemeyer ,  E.H.,  "Viscous  Corrections  to  Stewartson' s  Stability 
Criterion ,"  US  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving 
Ground,  Maryland,  Report  No.  1325 ,  June  1966.  (AD  489687 )  (See  Also  AG  ARP 
Conference  Proceedings ,  No.  10,  Mulhouse ,  France,  pp.  103-120 ,  September 
1966) 

4.  Murphy,  C.H.,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 
Liquid  Payload,"  US  Army  Ballistic  Research  Laboratory ,  Aberdeen  Proving- 
Ground,  Maryland,  ARBRL-MR-03194 ,  August  1982.  (AD  A118676)  (See  Also 
Journal  of  Guidance,  Control  and  Dynamics,  Vol.  6,  July-August  1983,  pp. 
280-286) 
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The  eigenvalue  Xg  exists  for  all  Re.  It  is  definitely  an  outlier  with 

respect  to  the  basic  set.  The  |Im  (x  )|  is  always  larger  than  that  for  any 

X^  and  its  inclusion  has  a  definite  effect  on  the  solution.  As  an 

example,  pressure  coefficient  at  the  intersection  of  the  endwall  and  sidewall 
is  41%  less  than  that  obtained  without  X  at  t  =  0.5  for  Re  =  10  and  A  = 
3.148.  s 

The  determination  of  X$  and  some  implications  of  its  inclusion  in  the 

theory  of  Reference  1  are  discussed  in  this  report.  The  simple  asymptotic 
methods  used  to  determine  Xg  show  their  efficacy  not  only  in  analysis  but  also 

in  the  successful  computation  of  the  pressure;  straightforward  calculation 
would  have  required  more  than  the  14  significant  figures  available  in  the  CDC 
7600. 


II.  PRECIS  OF  REFERENCE  1 

Certain  parts  of  the  theory  of  Reference  1  are  needed  here;  the  major 
points  are: 

a)  The  motion  of  the  cylinder  (projectile)  is  given  by  KQ  exp  (if$t) 

where  t  is  time,  f  =  (1  -  i <S )t ,  t  is  nutational  frequency/^  and  <5  is  the  yaw 
growth  rate.  KQ  is  the  magnitude  of  the  yaw,  assumed  small,  and  is  the 

parameter  used  to  linearize  the  Navier-Stokes  equations  in  the  perturbation 
analysis.  Lengths,  time,  velocity  and  perturbation  pressure  are  made  non- 

dimensional  by  a,  $_1,  a$  and  1'0  pa2$2,  respectively,  where  p  is  the  liquid 
density. 

b)  The  flow  variables  are  assumed  to  vary  as  exp  [i(f$t-9)]  where  9  is 
the  azimuthal  angle  coordinate;  the  variations  with  x  and  r,  the  longitudinal 
and  radial  coordinates,  respectively,  are  governed  by  linear  partial  differ¬ 
ential  equations  and  the  no-slip  conditions  on  the  sidewall  and  endwalls. 

c)  These  equations  are  solved  using  a  modal  analysis.  The  x- variation 
is  determined  by  a  second-order,  non-self-adjoint  system  for  spatial  eigen¬ 
values,  X,  and  eigenfunctions  sin  Xx  and  cos  Xx. 

d)  The  solution  satisfies  the  boundary  conditions  on  the  sidewall  but 
not  on  the  endwalls.  In  the  neighborhood  of  the  latter,  it  must  be  taken  as 
the  outer  solution  and  then  matched  to  an  inner  (boundary  layer  type)  solu¬ 
tion.  The  matching  is  asymptotic  for  Re  *  ». 

e)  First  term  matching  provides  a  solution  to  0  (Re-1)  and  corrected 
endwall  boundary  conditions  that  determine  the  X. 


f)  The  X  are  the  roots  of 


F  (z,e)  =  cos  z  +  ez  sin  z  =  0 


(2.1) 


where  z  =  A A,  e  =  6c/A  and  6c  is  the  complex  displacement  thickness  given  by 


6c  =  -(l+x)/2a(l-x)  +  (3-t)/23(1-t) 


a  =  [Re  (3-t )/2]i/^  (1- i ) 


(2.2) 


6  =  [Re(l+x)/2]*  (1+i ) ,  where  6  =  0,  f  =  t,  for 

convenience. 

The  range  of  x  is  limited  to  -1  <  t  <  3  in  order  that  inertial  waves  exist  in 
the  rotating  fluid;  this  is  the  range  in  the  inertial  coordinate  system 
employed  here  compared  to  the  frequency  range  (-2,2)  in  a  rotating  coordinate 
system. 

Since  (2.1)  is  solved  by  iteration,  first  guesses  for  the  roots  are 
required;  for  Re  >>  1 


z  -  (kn/2) (1-e)"  ,  k  odd, 


(2.3) 


was  used.  This  approximation  can  be  derived  in  several  ways;  see,  e.g.  (3.2). 

By  definition  t  is  real.  In  principle  there  is  no  restriction  on  t, 
except  t  *  +  1,  3.  In  practice  0  <  x  <  0.2  for  most  projectiles  and 
-0.5  <  t  <  0.5  for  most  experimental  simulations  of  a  spinning  projectile. 
Here  for  definiteness,  the  range  of  x  will  be  taken  to  be  -1  <  x  <  1.  From 
(2.2)  it  follows  that 


-xt /2  <  arg  e  <  -xr/4 


(2.4) 


for  -1  <  x  <  1. 


The  question  of  how  small  Re  can  be  and  still  obtain  reasonably  accurate 
solutions,  using  References  1  and  2,  will  not  be  discussed  because  it  is  not 
relevant  to  solving  (2.1). 


III.  ASYMPTOTIC  ANALYSIS  OF  THE  ROOTS 


A.  Limit  e  ♦  0 


The  roots  of  (2.1)  depend  on  the  single  parameter  e  which  is  proportional 

to  [ (1  - t jRe1  /4: ]~ 1 .  The  limit  of  primary  interest  is  e  +  0,  which  could 
correspond  to  Re  +  °°,  x  fixed.  The  limit  e  +  ”,  corresponding  to  x  +  1  and  Re 
fixed,  will  also  be  considered;  the  limit  x  -*■  -1  (or  even  x  -*■  3)  can  be 
treated  by  the  same  techniques  but  will  not  be  considered  here. 


>a 

.yi 


If  z  is  a  root  of  (2.1),  -z  is  also.  Since  the  eigenfunctions  are 
sin  Xx  and  cos  Xx,  and  because  of  the  symmetry  properties  of  the  solution1, 

only  the  roots  of  one  sign  are  required.  The  lesjlts  obtained  below  show  that 
only  roots  with 


-tt  /2  <  arg  z  <  0 


(3.1) 


are  appropriate.  First  let  z  =  0(1)  for  e  +  0  so  that  cos  z  is  the  dominant 
term  in  (2.1).  Thus  F  -*■  cos  z  and  z  -*■  +  kir/2,  k  =  1,3,...,  but  only  the 
positive  sign  need  be  considered.  This  gives  the  zeroth  approximation  to  the 
basic  set  of  roots, zj, ,  defined  by  zk  =  0(1)  or  zk  =  kir/2  for  e  ■*  0.  The  basic 

set  was  used  in  References  1-4.  The  zeroth  approximation  can  be  improved  by  a 
regular  perturbation.  Assuming  a  power  series  in  e  yields 


Z^  =  (kn /2) (1  +  e  +  . . . ) 


(3.2) 


so  that  -n/2  <  arg  z^  < 


It  is  not  possible  to  have  ez  sin  z  the  dominant  term  in  (2.1).  The 
possibility  that  neither  term  dominates  requires  ez  =  0(1).  Since  z  = 
0(l/e),  the  asymptotic  forms  for  sin  z  and  cos  z  are  required.  Using 


sin  z  =  (-i/2)(eiz  - 


=  (-i/2)eiz(l  - 


cos  z  =  (l/2)(eiz  +  e"iz) 


=  (l/2)e1z(l  +  e-2iz) 


(3.3) 


the  asymptotic  forms  are,  for  -n<  arg  z  <  0,  |z|  +  ®, 


sin  z  ~  ( - i / 2 ) e 1 


cos  z  ~  (l/2)e1z. 


(3.4) 


Substituting  (3.4)  into  (2.1)  yields  zQ,  the  zeroth  approximation  to  the 
asymptotic  solution  of  (2.1): 


zn  =  -i/e. 


(3.5) 


(For  0  <  arg  z  <  it,  zq  =  +i/e  is  obtained,  which  need  not  be  considered.) 

Note  that  because  of  (2.4),  -n/2  <  arg  zQ  <  0.  Therefore  zQ  =  O(Re^)  for 

Re  -►  00  and  is  distinct  from  any  member  of  the  basic  set.  Only  z  =  0(1)  and 

z  =  0(l/e)  roots  are  possible;  all  other  orders  give  contradictory  or 

inconsistent  results. 

The  eigenvalue,  which  is  O(Re^)  for  Re  00 ,  is  denoted  by  zs  -  z0. 

With  respect  to  the  basic  set  z$  is  a  singular  perturbation.  Depending  on  the 

iterative  method  used  to  solve  (2.1)  zs  can  be  determined  for  small  Re  (say 

Re  =  10)  and  x  >  .5  using  certain  0(1)  first  guesses.  However,  using  zQ  as 

the  first  guess  and  Newton's  method  to  solve  (2.1)  will  yield  zs  for  all  x  and 

Re  unless  Re  is  so  large  that  je12of  exceeds  the  capacity  of  the  computer  so 
that  the  Newton  method  fails.  But  for  such  Re,  z0  is  a  very  close 

approximation  to  zg.  Furthermore  the  approximation  zs  -  zQ  can  be  easily 

improved  so  that  the  iterative  solution  of  (2.1)  for  zs  is,  in  fact,  no  longer 
necessary  for  |e|  small  enough. 

Let  the  method  of  undetermined  coefficients  assume  a  series  expansion 

zs  =  zo(1  +  C1Y  +  (3-6) 

in  a  small  parameter  y.  With  (3.3),  (2.1)  can  be  written  as 

g(z)  =  1  +  e-2iz  -  iez(l  -  e'2iz)  =  0.  (3.7) 

Since  g  (zQ)  =  2e7^^ ,  the  small  parameter  y  =  e~3^c  is  suggested.  If  (3.6) 

contains  only  powers  of  y,  substitution  into  (3.7)  leads  to  a  contradiction  so 
that  the  ,  i  >  1,  cannot  be  determined.  The  reason  for  this  is  that  In  y 

terns  are  required  in  the  expansion.  This  can  be  shown  if  iteration  is  used 
to  solve  (3.7),  after  rewriting  it  as 


The  iteration  is  defined  by 


o 


l 


(3.8) 


■n+l 


;n  cn 

(1  +  Y  n)/(l  -  Y  n) 


and  the  result  for  n  =  1  is 


2s/z0  “  ?2  =  *  +  ^Y  +  2y2  +  ^Y2  1n  Y  +  0(y3  In  y), 


(3.9) 


The  method  of  undetermined  coefficients  works  with  these  functions  of  y. 
Equation  (3.9)  gives  the  first  four  terms  in  an  asymptotic  series  which 

appears  to  diverge  for  |y|>  10'^.  By  a  theorem  in  Reference  5  the  iteration 

(3.8)  converges  to  a  unique  root  if  |y|  <  10'  ,  approximately. 

Since  z$  and  the  basic  set,  z^,  are  unique  and  since  only  roots  of  0(1) 
and  0(l/e)  are  permissible,  all  roots  of  (2.1)  have  been  determined  for  e  0. 


B.  Limit  e  ♦  00 


Next  the  limit  e  »  is  considered,  i.e.,  t  -*•  1,  Re  +  0  or  both.  Because 
the  present  work  is  based  on  the  Re  -►  00  theory  of  Reference  1,  it  would  be 
i nappropri ate  to  allow  Re  0.  From  (2.2),  if  x  =  1-s  where  s  >  0  and  s  <<  1, 
then 


e  =  (-i/A  Re1/2s)[l  +  (3i s/4)  +  0(s2)].  (3.10) 


In  this  section  all  limits  are  taken  with  respect  to  e  -*■  00 .  Again  there  are 
only  two  kinds  of  roots  of  (2.1). 

(i)  z  =  0(1)  for  which  ez  sin  z  is  the  dominant  term  and  the  zeroth 
approximation  to  the  roots  is 


Zfj  =  Ntt  ,  N  =  1 ,  2 ,  . . . 


(3.11) 


5.  I saacson ,  E . ,  and  Keller,  H,  B.,  Analysis  of  numerical  Methods,  John  Wiley 
6  Sons,  New  York,  1966,  p.  86. 


J?~L  fJ. 


■'  ■  ■  (V  V. \V  i  vvr-- ' 
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the  root  z  =  0  being  discarded  since  it  is  not  0(1).  The  next  approximation 
gives 


zr,  =  Hw  -  (1/Nire) . 


(3.12 


This  expression  can  be  obtained  by  substituting  z  =  Nu  +  *  into  (2.1), 
expanding  sin  $  and  cos  $  in  power  series,  and  retaining  only  terms  through 
first  order  in  *. 

-1/2 

(ii)  z  =  0(c“  '  )  for  which  neither  term  in  (2.1)  is  dominant.  There  is 
a  single  root,  z^;  its  zeroth  approximation  is 


z'  =  -i/^/2 
so 


(3.13) 


with  -n/2  <  arg  ex/  <  0.  Although  it  has  not  been  proved  it  seems  clear  that 
zs  =  zs"  The  next  aPPr0xirnat10n  gives 


z'  =  (-i/e1/2)  -  ( i /3e3/2 ) . 


(3.14) 


The  approximations  (3.11)  and  (3.13)  or  (3.12)  and  (3.14)  can  be  used  as 
first  guesses  in  an  iterative  scheme  to  get  the  roots  of  (2.1);  for  the 
largest  e  tried,  |e|  =  1.08,  they  approximated  the  roots  to  only  one 

significant  figure.  Another  way  to  obtain  a  reliable  first  guess  for 

moderate  e  is  to  rewrite  (2.1)  in  terms  of  cot  z  (take  z  =  Nn  +  <t>)  and  use 
the  first  two  terms  in  the  expansion  of  cot  <j>.  The  result  is 


z  =  [-Nit  e  +  {(Nwe)2  -  4(e  -  1/3)  }1/2]/[2(e  -  1/3)]  +  Mtt  , 


(3.15) 


Some  representati ve  results  are  shown  in  Figures  1  and  2.  The  roots  of 
(2.1)  are  plotted  in  the  (zR,  Zj)  plane  for  Re  =  10  and  100,  x  =  0.1  and  A  = 

3.148  in  Figure  1.  All  the  zk,  k  =  1,  3,  ...  lie  close  to  the  real  axis; 

the  zs  for  the  two  values  of  Re  are  noted.  As  Re  increases,  the  z^  do  not 

change  much  but  zs  =  0(Re^2).  For  Re  =  10,  e  -  0.0836  -  0.1346  i  and 

Y  =  (-0.348  +  1 . 23 i )10-3  and  for  Re  =  100,  e  =  0.0264  -  0.0462  i  and 

Y  =  (-5.65  -  4.35i)10"^.  In  Figure  2,  the  roots  are  plotted  for  Re  =  10  and 

A  =  3.143  for  x  =  0.1  and  0.9.  For  the  latter  either  (3.5)  or  (3.13)  can  be 
used  as  first  guess  in  Newton's  method  to  obtain  zs  although  (3.13)  is  a 

better  guess,  as  expected.  Also  for  this  case  Real  (zs)  <  Real  (z^)  for  all 
odd  k. 
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IV.  SOME  EFFECTS  OF  INCLUDING  zs  IN  THE  THEORY  OF  REFERENCE  1 


The  existence  of  zs  has  a  number  of  consequences  in  the  theory  of 

Reference  1,  some  of  which  are  mentioned  here.  The  effect  on  the  non- 
dimensional  endwall  pressure,  Cp  at  r  =  1,  of  including  zs  is  shown  in  Table  1 

for  x  =  .5,  A  =  3.148,  and  four  values  of  Re.  Ten  terms  were  used  in  the 

series  to  calculate  Cp,  i.e.,  ten  values  of  z^  or  nine  plus  zs.  For  Re  <  10 

inclusion  of  zs  makes  a  significant  difference.  Since  the  theory  of  Reference 

1  is  asymptotic  for  Re  -*■  °°,  there  is  a  separate  question  of  how  accurate  that 

theory  is  for  Re  <  10°;  this  question  is  not  discussed  in  detail  here. 


Table  1.  Cp  at  r 


T  = 

.5 

A  =  3.148 

Re 

10 

102 

103 

10' 

With  z$ 

.0819 

.836 

.668 

.470 

Without  z,. 

.1398 

.793 

.662 

.470 

Figure  3  shows  Cp  vs  r  for  Re  =  100,  A  =  3.148  and  x  =  .1,  including  zs 

and  not.  The  effects  are  negligible  for  0  <  r  <  .45;  for  r  =  1  the  Cp  that 

includes  zs  is  16%  larger  than  that  if  zs  is  not  included.  The  results  for  Cp 

from  the  spatial  e.v.  method  (not  discussed  here)  are  also  shown;  it  can  be 
shown  that  these  are  accurate  to  the  number  of  significant  figures  plotted  and 
therefore  show  the  error  in  the  results  of  Reference  1,  even  including  z$. 

In  Reference  1,  the  satisfaction  of  the  sidewall  boundary  conditions, 
at  r  =  1,  requires  the  expansion  of  x,  -A  <  x  <  A,  in  terms  of  the 

eigenfunctions  sin  zkx/A  and  sin  zsx/A.  Two  methods  of  determining  the 

coefficients  were  given:  a  least  squares  fit  of  the  partial  sum  to  x  and  the 
eigenfunction  expansion  obtained  from  the  boundary  value  problem  satisfied  by 
the  e.v.  and  eigenfunctions.  This  is  not  a  self-adjoint  problem  so  the 

eigenfunctions  are  not  orthogonal  but  they  are  biorthogonal  with  respect  to 
the  solutions  of  the  adjoint  problem.  This  fact  enables  b^  and  bs,  the 

coefficients  in  the  expansion 


x  =  T  bk  sin  zk  x/A  +  bg  sin  z,.  x/A 


(4.1) 


mm® 


to  be  determined  analytically,  using 


bk  =  b  (zk) 

=  b(z$)  (4.2) 

b(z)  =  (2A/z2)  (1  +  e2z2)  sin  z/(l  +  e2z2  -  e). 

The  calculation  of  is  straightforward  because  z^  is  0(1)  but  the 

evaluation  of  bs  is  tedious,  time  consuming  and  could  require  more  than  the  14 

significant  figures  available  in  the  CDC  7600  for  Re  >  100.  The  reason  for 
this  is  that  the  numerator  cf  b$  goes  to  zero  exponentially,  which  can  be  seen 

from  (3.5),  (j.9)  and  the  fact  that  1  +  e2zs2  ~  -4y  =  -4e“2/,£.  If  bs  is 

evaluated  from  (4.2),  without  using  asymptotic  formulas,  a  large  number  of 
significant  figures  is  required  in  zs;  this  is  illustrated  in  Table  2  for  the 

quantity  called  u$(l)  =  -i  [2t  (1-t)/(1+t)]  b$  in  Reference  1. 

Table  2.  Effect  of  Significant  Figures  in  zs 


Re  =100  A  =  3.148  x  =  .1 


MSIG(zs) 

5 

7 

9 

11 

13 

Real (us(l) ) 

1.55  x  IQ-3 

1.63  x  IQ"5 

-2.85  x  io'6 

-2.12  x  10’6 

-2.132  x  io-6 

NSIG  (zs)  is  the  number  of  significant  figures  in  z$  used  in  (4.2);  similar 
behavior  is  found  for  Imag  (us  (1)).  In  contrast,  if  (3.5)  and  only  the  first 
two  terms  of  (3.9)  are  used  to  derive  the  asymptotic  form  of  bs, 

bs~  i  4A  e  e  '1/e, 

then  Real  (u  (1))  =  -2.131  x  10‘6  is  obtained  from  this  simple  calculation 

using  four  significant  figures  for  c  (i.e.,  z$);  this  differs  by  1  out  of 

2000  from  the  results  given  in  Table  2  for  NSIG  (zs)  =  13.  The  power  of 

asymptotic  methods  is  shown  by  this  example.  In  addition,  for  Re  >  100,  NSIG 
>  14  would  be  necessary,  at  which  point  it  would  be  impractical  to  use  the  C0C 
7600.  The  way  to  proceed  is  to  make  use  of  asymptotic  methods. 


A  check  on  the  relative  contribution  of  zs  can  be  made  by  computing  the 
sun  and  the  zs  terms  in  (4.1).  For  x  =  A  these  should  add  to  A.  With  the  sum 


denoted  by 
t  =  .1: 

£k,  tne  following  results  are 

obtained  for  Re  = 

<c 

o 

o 

r—i 

Real 

(£k)  =  2.9815 

Im 

<tk) 

=  .26793 

Real 

(b$  sin  z$)  =  .1665 

Im 

(bs  sin  zs) 

=  -.26799 

3.1480 

-6  x  10 

The  exact  results  are  3.148  and  0;  49  terms  were  used  in  Clearly  it  is 

essential  to  include  the  z$  terms. 

Since  bs  sin  z$  =  0  (e),  the  contribution  of  the  z$  term  to  the  series 

expansion  and,  ultimately,  the  pressure,  approaches  zero  as  Re  Since 

the  calculation  in  References  1  and  2  were  made  for  Re  >  1,000,  inclusion  of 
the  zs  term  would  have  a  negligible  effect  on  those  calculations. 
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LIST  OF  SYMBOLS 


radius  of  cylinder  [cm] 


c/a  =  aspect  ratio  of  cylinder 

coefficients  in  biorthogonal  eigenfunction  expansion  of 
f(x)  =  x  (See  (4.1)  and  (4.2)) 

half-height  of  cylinder  [cm] 

pressure  coefficient 

abbreviation  for  eigenvalue 

(1-i  5)  t 


cos  z  +  ez  sin  z 


odd  integer 


amplitude  of  projectile  angular  motion  [non-dimensional] 


radial  coordinate  /a 


Reynolds  number  =  a  $/v 


1  -  T 


time  [sec] 


■  i  [2f(l-T)/(l+T)]bt 


longitudinal  coordinate/a 


root  in  basic  set  of  roots  of  (2.1) 


root  in  basic  set  of  roots  of  (2.1)  for  e  +  ® 


LIST  OF  SYMBOLS  (Continued) 


0(e)  root  of  (2.1) 

0(e'1//z)  root  of  (2.1) 

-i/e1/2,  (see  (3.13)) 

-i/e,  (see  (3.5)) 
e-2A 

(1/t)  x  yaw  growth  per  radian  of  nutation 

complex  displacement  thickness 

Sc/A 

azimuthal  coordinate 
spatial  eigenvalue 

2 

kinematic  viscosity  of  liquid  [cm  /s] 
density  of  liquid  [g/cm3] 
nutational  frequency  of  cylinder/I 
spin  rate  of  cylinder  [radians/s] 
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